Molecular Plant

Sacha inchi (Plukenetia volubilis L.) is an emerging woody oilseed crop prized for its high alpha-linolenic acid (ALA) content. Despite its nutritional and economic value, the lack of high-quality genomic resources has hindered genetic improvement and the elucidation of its unique polyunsaturated fatty acid and lipid biosynthetic pathways. In this study, we report a high-quality, chromosome-scale genome assembly of sacha inchi with a total length of 710.62 Mb, integrated from Illumina, PacBio, and chromosome conformation capture (Hi-C) technology. The genome harbors 37,570 protein-coding genes, and 379.86 Mb (53.45%) of repetitive sequences. Phylogenomic analysis reveals that sacha inchi diverged from its closest relative Ricinus communis, [~] approximately 36.2 million years ago. Comparative genomics indicates that sacha inchi experienced only ancient whole genome duplication events. To elucidate the mechanisms governing ALA biosynthesis and triacylglycerol (TAG) accumulation in sacha inchi seeds, we performed temporal transcriptome profiling across six seed development stages. Our findings demonstrate that high TAG content is primarily driven by the sustained expression of biosynthetic genes and low activity of degradation genes during mid-to-late seed development. Notably, while genes encoding stearoyl-ACP desaturases (SADs) maintain the precursor pool, the expression of genes encoding fatty-acid desaturase 2 (FAD2) and fatty-acid desaturase 3 (FAD3) is positively correlated with the final accumulation of C18:2 and C18:3 fatty acids. We also identified lncRNAs as potential epigenetic regulators of these key pathways. This high-quality genome provides a critical foundation for elucidating the molecular mechanisms of seed growth and development in sacha inchi.

Plant amino acids function as both pathogen nutrients and essential drivers of systemic immunity. The regulation of amino acid homeostasis through transporters is a essential for mounting a robust and coordinated immune response in plants during pathogen infection. Using systems biology and integrative network science, we investigated bacterial virulence in Arabidopsis. By comparing gene coexpression networks of effector-triggered susceptibility (ETS) and pattern-triggered immunity (PTI), we uncovered a plant amino acid-related processes specifically linked to ETS. Integrating time-series transcriptomics, protein-DNA interactions, and mathematical simulations, we identified ANAC046 as a transcriptional regulator of amino acid processes during ETS. Single-cell RNA-Seq revealed that amino acid transporters are primarily expressed in companion and mesophyll cells, while functional validation confirmed ANAC046s roles in promoting susceptibility. Further integration of transcriptome and interactome data showed that amino acid-related genes interact with key immune hub proteins. Network topology analysis enabled the characterization of seven additional genes involved in plant defense. To support community-wide research, we developed MIData, an open-access platform for pre-analyzed Arabidopsis networks. Together, our findings demonstrate the power of systems-level approaches in uncovering hierarchical regulatory mechanisms underlying plant susceptibility to bacterial pathogens.

Hi-C data is commonly used for reference-free de novo scaffolding. However, with the rapid increase in high-quality reference genomes, reference-guided workflows are now more practical for assembling large numbers of target genomes without relying on costly and labor-intensive Hi-C sequencing. Recently, a pangenome graph-based haplotype sampling algorithm was introduced to generate personalized graphs for target genomes. Such graphs have strong potential as references for reference-guided contig scaffolding. Here, we present noHiC, a reference-guided scaffolding pipeline supporting key steps of plant contig scaffolding. A distinctive feature of noHiC is the nohic-refpick script, generating a best-fit synthetic reference (synref) from a pangenome graph that is genetically close to the target contigs. This enables the integration of genetic information from many references (up to 48 in our tests) without using them separately during scaffolding. Synrefs showed advantages over highly contiguous conventional references in reducing false contig breaking during reference-based correction. Additionally, nohic-refpick can be combined with fast scaffolders (ntJoin) to rapidly produce highly contiguous assemblies using synrefs derived from pangenome graphs. The noHiC pipeline, used alone or in combination with ntJoin, can generally produce assemblies that are structurally consistent with public Hi-C-based or manually curated genomes. The pipeline is publicly available at https://github.com/andyngh/noHiC. GRAPHICAL ABSTRACT O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=82 SRC="FIGDIR/small/712436v1_ufig1.gif" ALT="Figure 1"> View larger version (9K): org.highwire.dtl.DTLVardef@40bd8forg.highwire.dtl.DTLVardef@5d2bbborg.highwire.dtl.DTLVardef@e214a3org.highwire.dtl.DTLVardef@b90b06_HPS_FORMAT_FIGEXP M_FIG C_FIG

Flexible developmental programs enable plants to customize their organ size and cellular composition. In leaves of eudicots, the stomatal lineage produces two essential cell types, stomata and pavement cells, and plants can adjust the total numbers and ratios of these cell types in response to external cues. Central to this flexibility is the stomatal lineage-initiating transcription factor, SPEECHLESS (SPCH). Here we explore the mechanisms underlying SPCHs involvement in environmental response. Using multiplexed CRISPR/Cas9 editing of SlSPCH cis-regulatory sequences in tomato, we identified variants with altered stomatal development responses to drought, light and temperature cues. By creating and live-cell tracking translational reporters of SlSPCH and its paralogues SlMUTE and SlFAMA, we revealed the corresponding cellular events that lead to the environmental change-driven responses in stomatal production and leaf form. Plants bearing the novel reporters and SlSPCH variants are powerful resources for fundamental and applied studies of tomato resilience in response to climate change.

Sesame (Sesamum indicum) is one of the earliest domesticated oilseed crops and is valued for antioxidant lignans that stabilize oil quality. However, the genomic and evolutionary history of the genus Sesamum, including the origin of its allotetraploid relative S. radiatum and the diversification of lignan metabolism, remains poorly understood owing to limited chromosome-scale genomic resources. Here we present chromosome-level genome assemblies for three wild Sesamum species, two Ceratotheca species and a Japanese sesame cultivar to reconstruct genome and karyotype evolution across the Sesamum-Ceratotheca complex. Comparative analyses show that the derived x=16 lineage originated from an ancestral x=13 karyotype through chromosome fission, fusion and translocation, whereas another x=13 lineage underwent extensive restructuring associated with retrotransposon expansion. Phylogenomics places Ceratotheca within the x=16 Sesamum clade and reveals that S. radiatum originated through hybridization involving a C. sesamoides-like ancestor. The antioxidative lignan gene CYP92B14 was reintroduced via the BB progenitor, linking hybridization with restoration of oil-stabilizing metabolism during sesame evolution.

The establishment of aphid-plant interaction involves the secretion of a salivary MIF protein. Morphological analyses revealed that aphid MpMIF1 prevents plant cell death, protects organelles from stress, and may promote plant cellular recovery. Co-expression of aphid MpMIF1 and the cell death inducer Npp1 revealed that MpMIF1 modulates autophagy-related genes ATG7/BECLIN1, impair plant senescence regulator ATAF1 and regulate apoptosis-like via Caspase-3-like activity. This effect on multiple-cell death pathways helps to maintain cellular homeostasis during aphid infection. Investigations on DNA Damage Response (DDR) signaling pathways demonstrated that aphid MpMIF1 reduces {gamma}H2A.X phosphorylation, maintains activity of the DNA repair protein RAD51 and stabilizes cell cycle checkpoint expression WEE1 under genotoxic stress. Therefore, MpMIF1 actively participates to the maintenance of a functional DDR. Finally, we showed that aphid MpMIF1 physically interacts with SOG1, a functional analog of animal p53 and central regulator of DDR, cell cycle arrest and programmed cell death in plants. These findings establish MpMIF1 as a key regulator of plant cell death during aphid-plant interactions and highlight its potential as a biotechnological tool for protecting major crops against aphid infection.

Strigolactones (SLs) are plant hormones regulating shoot branching and symbiotic interactions, but their trace-level abundance limits research and applications. Here, we optimized a Nicotiana benthamiana transient expression system for SL production by tuning agroinfiltration parameters and co-expressing rate-limiting carotenoid biosynthetic genes. Overexpression of Zea mays PSY1 or an Arabidopsis PSY-GGPS11 fusion increased carlactone production over 2-fold and enhanced downstream SL accumulation. Using this platform, we discovered that sorghum cytochrome P450 SbCYP728B35 catalyzes conversion of 5-deoxystrigol to sorgolactone, revealing a previously unknown function. These results establish metabolic engineering of precursor supply as an effective strategy for boosting SL production and demonstrate N. benthamiana as a robust system for pathway elucidation and biotechnological synthesis of bioactive strigolactones.

Salt stress alters plant development, including the floral transition, but regulation of timing of flowering by salt is poorly understood at the molecular level. To identify genetic loci regulating the floral transition under high soil salinity, we performed a genome-wide association study (GWAS) in Arabidopsis thaliana and identified natural variation at the UGT74E1-UGT74E2-BT3 (UUB) locus that correlates with bolting time specifically in response to salt stress. Genetic analysis revealed BT3 as a novel repressor of the floral transition in control conditions. Similarly, the putative IBA glycosylases UGT74E1 & UGT74E2 delay the floral transition in control conditions. Furthermore, we identified that IBA homeostasis regulators TOB1 and ECH2/IBR10 play a key role in the floral transition, and that ECH2/IBR10 are required for the early flowering phenotype of the ugt74e1/ugt74e2 double mutant, indicating that UGT74E1 & UGT74E2 delay flowering by altering IBA homeostasis. A pangenome analysis of the UUB locus revealed variation in the occurrence of the DNA transposon SAUERKRAUT (SKRT). CRISPR-mediated SKRT deletion in Col-0 affected gene expression both within and outside the UUB locus and caused a salt-dependent delayed floral transition. The delayed bolting phenotype of the skrt-2 mutant also depends on ECH2/IBR10 function, indicating that SKRT accelerates the floral transition by altering IBA homeostasis. Finally, targeted demethylation of SKRT resulted in delayed floral transition under salt stress. Taken together, our data show a role for SKRT and its DNA methylation levels in the salt-dependent bolting time response in Arabidopsis, revealing a novel molecular mechanism to control flowering in adverse conditions.

Concerns over climate change have intensified the demand for stress resistant crops like hybrid wheatgrass (HWG; Elymus hoffmannii, StStStStHH), a perennial forage species known for its exceptional salt and drought tolerance. However, hexaploidy and high heterozygosity have complicated efforts to resolve its genomic structure and evolutionary history. Here, we present high-quality, haplotype-resolved, chromosome-level genome assemblies for HWG (CDC Saltking) and its putative progenitor, bluebunch wheatgrass (Pseudoroegneria spicata). By integrating PacBio HiFi and ultra-long Oxford Nanopore sequencing with Hi-C scaffolding, we assembled the 10.7 Gb HWG genome into 21 pseudochromosomes per haplotype. Our phylogenomic analysis redefines the origin of the H subgenome, positioning it as an intermediate between Old-World Hordeum marinum (sea barley) and Hordeum brevisubulatum. Notably, we identified significant chromosomal rearrangements, including a unique duplication on St chromosome 4. Transcriptome analysis across multiple tissues revealed a pronounced expression dominance of the H subgenome. This dominance was not associated with reduced LTR density, suggesting that selective pressures for rapid adaptation of the latest subgenome entrant may drive its dominance. Finally, using the f-branch statistic, population genomic analysis of 189 accessions representing eight Elymus and Pseudoroegneria species revealed extensive reticulate evolutionary relationships and identified P. spicata as a major, asymmetric genetic donor within the wheatgrass complex. These resources provide a foundational framework for future genomic research and genetic improvement in grasses and for the introgression of stress-tolerance traits into cereal crops such as wheat. Key MessagesDevelopment of world-first high-quality chromosomal-level haplotype-resolved genome assemblies of hexaploid HWG and diploid progenitor, Pseudoroegneria spicata, enabled the identification of the subgenome origins. This study resolved the evolutionary placement of the St genome and clarified the history of polyploidization and hybridization in HWG. Homeolog expression bias in the H subgenome likely reflects selective pressure favoring greater gene retention and upregulation of functionally important genes, thereby enhancing hybrid fitness. Population structure analysis distinctly differentiates P. spicata, E. repens, E. hoffmannii from other European Pseudoroegneria species. The findings reveal the complex patterns of interspecific gene flow and population dynamics within the Elymus and Pseudoroegneria species.

Durable resistance to soil-borne pathogens remains elusive in cereals, partly because susceptibility (S) genes that facilitate root infection have not been identified in monocots. In the model legume Medicago truncatula, the SCAR/WAVE complex member MtAPI functions as a root S-gene for microbial invasion. Whether SCAR gene associated susceptibility function is conserved in monocots, and whether SCAR gene inactivation can enhance root resistance in cereals, remains unknown. Here, we identify and characterize three SCAR genes in barley: HvSCAR-A, HvSCAR-B, and HvSCAR-C. Cross-species complementation assays indicate that HvSCAR-B and HvSCAR-C are functionally similar to MtAPI. While hscar-b and hvscar-c single mutants exhibited no major growth defects, hvscar-a mutants showed strongly reduced seed production, and a hvscar-b/c double mutant displayed shorter root hairs. Notably, the hvscar-b/c double mutant exhibited increased resistance to the hemibiotrophic pathogen Phytophthora palmivora but greater colonization by the symbiotic arbuscular mycorrhizal fungus Funneliformis mosseae, underscoring a complex role in plant root - microbe interactions. Our findings reveal a conserved susceptibility function of SCAR genes in monocots and identify api monocot homologs as promising targets for engineering disease resistance in cereals. This study offers new insights into SCAR protein functional diversification and its potential for improving root health in crop plants.

Pepper (Capsicum annuum L.) is a recalcitrant species regarding shoot regeneration, a trait that serves as a major bottleneck for the application of genetic engineering tools. In this study, comparative genetic analysis between a rare high-regeneration cultivar and a common low-efficiency cultivar identified a single nucleotide polymorphism (SNP) in PHYTOCHROME A SIGNAL TRANSDUCTION 1 (CaPAT1) that determines shoot regeneration efficiency. The T478C SNP in the high-efficiency cultivar converts a stop codon into an Arg codon, leading to translational read-through into the neighboring gene and forming an intact GRAS domain. This SNP-mediated formation of full-length CaPAT1 is essential for its dimerization. Notably, the overexpression of CaPAT1T478C in multiple low-efficiency cultivars, including both hot and bell peppers, significantly improved both shoot regeneration and transformation efficiency in the transformed T0 generation. These findings demonstrate the pivotal role of CaPAT1 in enhancing shoot regeneration and provide a robust strategy to overcome recalcitrance in pepper.

Macroautophagy is a conserved intracellular catabolic process in eukaryotes that participates in chloroplast degradation, through the selective breakdown of chloroplast components. Selective autophagy of membrane-bound organelles typically requires receptors that bridge organelle membranes and pre-autophagosomal structures. Here we identify OEP24.1 as a new receptor in the selective chloroplast piecemeal autophagy, supporting the degradation of stromal proteins. We found that the {beta}-barrel protein OEP24.1 is located at the outer membrane of plastid envelopes and on bodies budding off plastids into the cytosol and containing stroma proteins. OEP24.1 interacts physically with ATG8 autophagy proteins in a UIM dependent manner. OEP24.1-GFP and RFP-ATG8 colocalize with in mobile autophagosome-like puncta in the cytosol and in autophagic bodies within the vacuole. Delivery of OEP24.1 to vacuole lumen is dependent on active autophagy. OEP24.1 controls carbon allocation at the whole plant level, carbon concentrations in flowering stems and xylem composition. These phenotypes can be explained by the role of OEP24.1 in metabolite diffusion across the chloroplast envelope, and by its involvement in the facilitation of chloroplast quality control through piecemeal autophagy.

Plant-type (Form I) Ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) suffers from inherent catalytic trade-offs and a strong dependency on other proteins--including an essential small subunit (SSU) and auxiliary chaperones--for assembly, constraining the enzymes evolutionary and engineering potential. Here, we investigated representatives from the newly discovered clade Form IF. These enzymes do not require specific chaperones to form functional complexes, exhibit high CO2-specificities (SC/O [~]50) while maintaining high turnover rates (up to kcat [~]11 s-1). Remarkably, two Form IF representatives (IF-1/IF-2) lost the dependency on the SSU and assemble into homo-octameric complexes without their cognate SSUs. While the SSU is not necessary for catalysis, its addition improves both activity and specificity in IF-1/IF-2. Our results show that complexity is actually not required to achieve highly active, specific and functional Rubisco variants--and that this complexity can even be reverted--which challenges our current thinking on the evolution and catalytic mechanism of Rubisco.

The global antibiotic resistome remains largely unexplored, not because antibiotic resistance genes (ARGs) are rare in the environment, but because many are evolutionarily distant from known ARGs. Current computational approaches primarily rely on sequence homology, and thus miss distant homologues. We develop GeoARG, a geometry-enhanced framework that integrates structural features with protein language models through knowledge distillation, enabling efficient large-scale screening using sequence input alone. Across multiple benchmarks, GeoARG substantially improves the detection of remotely homologous ARGs, particularly under low sequence identity and fragmented conditions. Large-scale metagenomic analysis uncovers 1,485 high-confidence ARG candidates that are highly divergent from known ARGs, expanding the phylogenetic and functional landscape of the resistome. Structural analyses further show that these candidates preserve active-site geometry and maintain stable ligand-binding configurations consistent with known resistance mechanisms. These results demonstrate that geometric constraints enable systematic expansion of the resistome and facilitate the discovery of evolutionarily distant yet functionally conserved genes. A public web server is available at https://ycclab.cuhk.edu.cn/GeoARG/.

The transition from hulled to free-threshing grain was a pivotal event in wheat domestication, enabling efficient harvesting and processing. Threshability in tetraploid wheat is controlled primarily by the Q locus and two Tenacious glume (Tg) loci on chromosomes 2A and 2B, yet the molecular basis of the major Tg1-B locus remains incompletely characterized. Here, we phenotyped a durum wheat x wild emmer wheat (WEW) recombinant inbred line (RIL) population across two field environments and performed QTL analysis for glume tenacity (TG), threshability ratio (THRR), and seed number per spike (SDNPS). A total of 19 significant QTLs were detected across six chromosomes. The largest-effect loci for both TG and THRR co-localized on chromosome 2B, with LOD scores up to 14.22 and phenotypic variance explained up to 31.2%, corresponding to the previously described Tg1-B locus. To validate this QTL, the donor RIL was backcrossed three times to Svevo to generate a near-isogenic line, NIL-65 (BC3F5), confirmed by whole-genome skim sequencing to carry a homozygous WEW introgression at Tg1-B. A segregating BC4F2 population derived from NIL-65 confirmed that plants homozygous for the dominant Tg1-B allele displayed significantly higher glume tenacity and intact glume morphology compared to tg1-B sister lines, which exhibited basal glume cracking characteristic of the free-threshing phenotype. Genotyping-by-sequencing delimited the causal interval to an approximately 11 Mb introgression on chromosome 2B. These results confirm the major role of Tg1-B in determining glume tenacity in tetraploid wheat, provide a validated near-isogenic germplasm resource, and lay the foundation for fine-mapping and functional characterization of the underlying gene(s).

Hundreds of proteins in the cell require iron (Fe) or Fe-containing cofactors to function. However, how Fe2+ or Fe3+ are specifically allocated to each of these proteins in plant cells remains largely unknown. It has been proposed that Fe metalation could be driven by specific interactions with Fe-shuttling proteins known as Fe-chaperones. Here, we present the first family of plant Fe2+-chaperones (ICHAPs) with orthologues in dicots and monocots. The role of these proteins in Fe distribution to Fe-dependent metabolic processes has been illustrated using symbiotic nitrogen fixation in Medicago truncatula root nodules. ICHAP1 is a soluble Fe2+-binding protein that interacts with plasma membrane Fe2+ transporter NRAMP1, but not with symbiosome Fe2+-transporters. ICHAP1 mutants present altered Fe distribution in cells and they cannot fix nitrogen. A second family member, ICHAP2 is required to target Fe2+ to symbiosomes, as it accepts Fe2+ from ICHAP1 and interacts with symbiosome Fe2+-importer VTL8, but not with NRAMP1. These results indicate a path for Fe2+ allocation from the plasma membrane to the symbiosome through specific protein-protein interactions and Fe2+ exchange from NRAMP1 to ICHAP1, to ICHAP2, and to VTL8.

The combination of target capture sequencing (TCS) with low-coverage whole genome sequencing (lcWGS), an approach known as Hyb-Seq, has allowed the integration of natural history collections into the genomics revolution, transforming biodiversity research. To implement Hyb-Seq, a collection of genomic targets, often nuclear orthologs, is needed to design probes for TCS. In flowering plants, the universal Angiosperms353 probe set has been proven resolutive at multiple evolutionary scales, with caveats. Malpighiales is known to be one of the most challenging flowering plant orders to resolve. Within this order, the clusioid clade ([~]2.2K species, 94 genera, five families) is no exception. To resolve phylogenetic relationships in this recalcitrant clade, we design a custom probe set, the Clusioids626 kit, composed of 39,936 120-mer probes targeting 626 nuclear orthologs ([~]6.6M nucleotides). This probe set includes all Angiosperms353 targets and 273 clusioid-specific ones, carefully chosen taking copy-number, length evenness, and phylo-informativeness into account. We test our probe set on 70 accessions representing all families and tribes in the clusioid clade. On average, 50.4% of TCS reads mapped to our targets, recovering a median of [~]600 orthologs. Relationships for all clusioid families are fully resolved for our nuclear targets. Additionally, 105 plastid coding DNA sequences were retrieved from the lcWGS fraction. A strong cyto-nuclear conflict was detected. The Clusioids626 kit performs better than the universal Angiosperms353 enrichment panel alone. Our kit design workflow can be extended into other lineages for which a universal probe set exists but more resolution is needed.

T-DNA insertion mutants are widely used to disrupt genes and infer their functions, yet the insertions can also trigger unintended genomic changes that confound phenotypic interpretation. Here, we used T-DNA insertion mutants affecting the major mitochondrial malate dehydrogenase (MDH1) and the heterodimeric NAD-dependent malic enzymes (ME1 and ME2) to examine their functional coordination across photoperiods and irradiance regimes. Under short days, especially at low light intensity, mdh1xme2 mutants were markedly smaller than wild type and, unexpectedly, than the mdh1xme1xme2 triple mutant, and they showed a more pronounced reduction in photosynthetic capacity. ME1 was undetectable in mdh1xme2, implying that the double and triple mutants effectively lack heterodimeric ME and should therefore behave similarly, contrary to what we observed. Whole-genome analysis resolved this discrepancy by revealing that the MDH1 T-DNA insertion in mdh1xme2 is accompanied by a major rearrangement, a 137-kbp duplication downstream of the insertion site, which was absent in the mdh1xme1xme2 triple mutant. This duplication increased gene dosage and elevated transcript abundance across the duplicated interval, while proteomics detected 5 of the 38 encoded proteins, including PEPC1. mdh1xme2 accumulated oxaloacetate-derived amino acids and displayed an altered carbon/nitrogen balance, making PEPC1 a plausible contributor to the exacerbated mdh1xme2 phenotype. Together, our data indicate that a T-DNA-linked structural variant can amplify expression of dozens of genes and intensify phenotypes at specific conditions, thereby affecting the interpretation of genotype-phenotype relationships. Because Agrobacterium-mediated DNA transfer also underpins many genome-editing workflows, our findings argue that structural validation around insertion/editing loci should be considered essential when interpreting T-DNA-derived plant lines.

Poplar seed fibers cause environmental and health concerns, yet their developmental mechanisms remain poorly understood. Here, we constructed a high-resolution spatiotemporal transcriptomic atlas of female poplar capsules by integrating single-nucleus and spatial transcriptomics. We delineated the developmental trajectory of seed fibers, confirming their origin from placenta cells, and identified three functionally distinct fiber cell subtypes involved in initiation, metabolic support, and elongation. Weighted gene co-expression network analysis (WGCNA) identified several hub transcription factors, including PtoMYB, PtoHDT1, PtoEIF6 and PtoPDF2, that may serve as key regulators of fiber development. Our study provides a cellular-resolution framework for understanding trichome development in woody perennials and offers candidate targets for functional characterization toward breeding low-fluff poplar cultivars. HighlightsO_LIA spatiotemporal transcriptomic atlas of poplar capsule development is constructed at single-cell resolution C_LIO_LIFiber cells originate from placenta cells and comprise three functionally distinct subtypes C_LIO_LIProvides molecular targets for breeding low-fluff poplar cultivars to mitigate environmental pollution C_LI

Plant mitochondrial genomes are evolutionary paradoxes, structurally fluid yet functionally conserved. In hybridization-prone lineages like Camellia (Theaceae), where polyploidy and introgression blur species boundaries, mitochondrial dynamics remain unexplored. Here, we assemble mitochondrial genomes of eight Camellia species and Actinidia eriantha as outgroup, revealing three hallmarks of reticulate evolution: first, genome inflation (788,876 bp-1,000,198 bp) driven by short repeats (<100 bp; [~]92.6-94.9% of total repeats, PGLS: R{superscript 2}=0.9729, P=6.293x10-6), overturning the long-repeat paradigm in plants; second, multichromosomal architectures in C. sinensis var. sinensis LJ43 and C. taliensis, with lineage-specific plastid DNA occupying up to 25.5% of mitochondrial content; and third, phylogenomic conflicts between mitochondrial, chloroplast, and nuclear gene trees, revealing genome-wide hybridization signatures. Mitochondrial genomes emerge as recorders of historical gene flow, capturing repeat explosions triggered by pervasive hybridization and cytonuclear co-evolution. By linking micro-repeat dominance to genomic shock, this work redefines mechanisms of plant mitochondrial evolution and provides a roadmap to resolve taxonomic complexities in rapidly radiating angiosperms.